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Abstract 

Although the physical properties of the 2D and ID Ising models are quite different, we 
point out an interesting connection between their complex-temperature phase diagrams. We 
carry out an exact determination of the complex-temperature phase diagram for the ID 
Ising model for arbitrary spin s and show that in the Us = e~^^'^^ plane (i) it consists of 
Nc,iD = 4,s^ infinite regions separated by an equal number of boundary curves where the 
free energy is non-analytic; (ii) these curves extend from the origin to complex infinity, and 
in both limits are oriented along the angles 6'„ = (1 + 2n)7r/(4s^), for n = 0, 4s^ — 1; (iii) 
of these curves, there are Nc^^e,id = ^c.atvk,!!) = [-^^l the first and second (NE and NW) 
quadrants; and (iv) there is a boundary curve (line) along the negative real Ug axis if and 
only if s is half-integral. We note a close relation between these results and the number of 
arcs of zeros protruding into the FM phase in our recent calculation of partition function 
zeros for the 2D spin s Ising model. 
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Recently, we presented calculations of the complex-temperature (CT) zeros of the parti- 
tion functions for square-lattice Ising models with several higher values of spin, s = 1, 3/2, 
and 2 In the thermodynamic limit, these zeros merge to form curves across which the free 
energy is non-analytic, and thus calculations for reasonably large finite lattices give insight 
into the complex-temperature phase diagrams of these models. These phase diagrams consist 
of the complex-temperature extensions of the Z2-symmetric, paramagnetic (PM) phase; of 
the two phases in which the Z2 symmetry is spontaneously broken with long-range ferromag- 
netic (FM) and antiferromagnetic (AFM) order; and, in addition, certain phases which have 
no overlap with any physical phase (denoted "O" for other). Some of the zeros lie along 
curves which, in the thermodynamic limit, separate the various phases. In addition, there 
are zeros lying along various curves or arcs which terminate in the interiors of the FM and 
AFM phase. Physical and CT singularities of the magnetization, susceptibility, and specific 
heat obtained from analysis of low-temperature series have been discussed recently for the 
square lattice Ising model with the higher spin values s = 1 and s = 1, 3/2, 2, 5/2, and 
3 1- 

Here we give some further insight into the complex-temperature phase diagrams for 
higher-spin Ising models. We first report an exact determination of the CT phase diagrams 
of the ID Ising model for arbitrary spin s. We then point out a very interesting connection 
between features of these ID phase diagrams and certain properties of the phase diagrams 
of the Ising model on the square lattice inferred from our calculation of partition function 
zeros for higher spin values. This connection is useful because, unlike the 2D spin 1/2 case, 
no exact closed-form solution has ever been found for the 2D Ising model with spin s > 1, 
and hence further elucidation of its properties is of continuing value, especially insofar as 
these constrain conjectures for such a solution. Of course, the physical properties of a spin 
model at its lower critical dimensionality (here (i^.c.d. = 1) are quite different from those for 
d > di c.d.- However, as we shall discuss, some of the properties of the CT phase diagram for 
d = 2 exhibit simple relations with the d = 1 case. [| 

There are several reasons why CT properties of statistical mechanical models are of 
interest. First, one can understand more deeply the behavior of various thermodynamic 
quantities by seeing how they behave as analytic functions of complex temperature. Second, 
one can see how the physical phases of a given model generalize to regions in appropriate 
CT variables. Third, a knowledge of the CT singularities of quantities which have not been 
calculated exactly helps in the search for exact expressions for these quantities. Fourth, one 
can see how CT singularities in functions such as the magnetization and susceptibility are 

^Indeed, from the d = 1 + e and d = 2 + e expansions for the Ising and 0{N) models [Q, one knows that 
expansions above di,c.d. can even give useful information about physical critical behavior. 
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associated with the boundaries of the phases and with other points where the free energy is 
non-analytic. Such CT properties were first considered (for the 2D, s = 1/2 square-lattice 
Ising model) in Ref. |^ and for higher-spin (2D and 3D) Ising models in Ref. 0. 

The spin s (nearest-neighbor) Ising model is defined, for temperature T and external 
magnetic field H, by the partition function Z = I]{5„}e~^^ where, in a commonly used 
normalization, 

n=-{J/s^) J2 SnSn'-{H/s)J2Sn (1) 

<nn'> n 

where Sn G {—s, —s + 1, s — 1, s} and (3 = (ksT)^^. H = unless otherwise indicated. We 
define K = jSJ and then a generalized (i.e. with negative as well as positive 

powers) polynomial in Ug- The (reduced) free energy is / = —(3F = limi^i^^oo In Z in 
the thermodynamic limit. 

For d = 1, one can solve this model by transfer matrix methods. One has 

2s+l 

Z = Tr{T^) = Y: (2) 

where the Xsj, j = 1, ...2s + 1 denote the eigenvalues of the transfer matrix T defined by 
Tnn' =< 'n\exp{{K/s'^)SnSn')\n' > (we assume periodic boundary conditions for definite- 
ness). It is convenient to analyze the phase diagram in the Us plane. For physical temper- 
ature, phase transitions are associated with degeneracy of leading eigenvalues 0. There is 
an obvious generalization of this to the case of complex temperature: in a given region of 
Us, the eigenvalue of T which has maximal magnitude, Xmax, gives the dominant contribu- 
tion to Z and hence, in the thermodynamic limit, / receives a contribution only from Xmax'- 
f = In(Amax)- For complex K, f is, in general, also complex. The CT phase boundaries are 
determined by the degeneracy, in magnitude, of leading eigenvalues of T. As will be evident 
in our ID case, as one moves from a region with one dominant eigenvalue Xmax to a region 
in which a different eigenvalue A[„„^ dominates, there is a non-analyticity in / as it switches 
from / = \n{Xmax) to / = \n{X'^^^). The boundaries of these regions are defined by the 
degeneracy condition |Ama^| = |A^„^|. These form curves in the Ug plane.0 

Of course, a ID spin model with finite-range interactions has no non-analyticities for any 
(finite) value of K, so that, in particular, the ID spin s Ising model is analytic along the 
positive real Us axis. For a bipartite lattice, Z and / are invariant under K —>■ —K, i.e.. 
Us — > 1/ms. It follows that the CT phase diagram also has this symmetry, i.e., is invariant 
under inversion about the unit circle in the Ug plane. This symmetry also holds for a finite 
bipartite lattice; for c? = 1, the lattice is bipartite iff is even, and for our comments 
^By "curves" we include also the special case of a line segment. 
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about finite-lattice results, we thus make this restriction. Further, since the Xgj are analytic 
functions of Ug, whence \sj{u*) = \sj{us)*, it follows that the solutions to the degeneracy 
equations defining the boundaries between different phases, lA^jl = \Xs,e\, are invariant under 
i^s ^ ^s- Hence, the complex-temperature phase diagram is invariant under Ug u-l. 

Wc shall present results for a few s values explicitly. For s = 1/2, one has {ui/2Y^^Xi/2,j = 
l±('Ui/2)^^^- / is an analytic function of Ui/2 except at points which constitute the solution to 
1^1/2, i| = 1^1/2,2!; these comprise the negative real axis, —00 < Ui/2 < 0. Apart from this line, 
the dominant eigenvalue of T is Ai/2,i- For s = 1, the eigenvalues of are Ai_i = u^^ — Ui 
and 

Ai,,=2,3 = (1/2) [th' + 1 + Ml ± {uf - 2u^^ + 11 - 2ui + ulf'^] (3) 

As shown in Fig. 1(a), the phase diagram consists of four phases, the complex-temperature 
extension of the PM phase, together with three O phases. The curves separating these 
phases are the solutions of |Ai,i| = |Ai.2|. The third eigenvalue, Ai^s, is always subdominant. 
In the two phases containing the real axis, Ai 2 has maximal magnitude, while in the 
two containing the imaginary Ug axis, Au is dominant. The CT zeros of Z calculated for 
finite lattices lie on or close to these curves, starting a finite distance from the origin and 
being distributed in a manner symmetric under the inversion Ug l/ug. As the lattice size 
increases, the zeros spread out, the one with smallest (largest) magnitude moving closer to 
(farther from) the origin. For s — 3/2, the eigenvalues are given by 

C[(l - 2u^ + Au^ + u^ + u^ + Au"" - 2u^ + u^) - 2r]u^/^{l - 6u^ + u")] (4) 

where here u = U3/2 and {rj, Q = {+, +), (-, +), (+, -), (-, -) for j = 1, 2, 3, 4. The phase 
diagram is shown in Fig. 1(b) and consists of nine regions separated by the curves where 
1-^3/2,11 = |A3/2,2|- In the region containing the positive real Ug axis, A3/2,i is dominant, and 
as one makes a circle around the origin, each of the nine times that one crosses a boundary, 
there is an alternation between A3/2,i and A3/2,2 as the dominant eigenvalue. 

We find the following results for general s: (i) the complex-temperature phase diagram 
consists of 

Nc,iD = (5) 

(infinite) regions separated by an equal number of boundary curves where the free energy is 
non-analytic; (ii) these curves extend from the origin to complex infinity, and in both limits 
are oriented along the angles 

(l + 2n)7r 
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for n = 0, ...4s^ — 1; (iii) of these curves, there are 

Nc,NE,lD = ^c,Ariy,lD = [s^] (7) 

in the first and second (NE and NW) quadrants, where [x] denotes the integral part of x; 
and (iv) there is a boundary curve (which in this case is a straight hne) along the negative 
real Ug axis if and only if s is half-integral. {Nc^se,id = Nc^ne,id and Nc^sw,id = -^c.Afw.iz? 
by the u u* symmetry.) To derive these results, we carry out a Taylor series expansion of 
the A^j in the vicinity of Ug = 0. The dominant eigenvalues have the form 



uf Xg = 1 + ... + ttg jui^^ + ... 



where the first ... dots denote terms which are independent of j and the second ... dots 
represent higher order terms which are dependent upon j. Setting Ug = re*^ and solving the 
equation \Xgj \ = \Xs,e\ yields cos(2s^6') = 0, whence 

2s'^e = ^ + mr, n = 0,... 45^-1 (9) 

Each of these solutions yields a curve across which / is non-analytic, corresponding to a 
switching of dominant eigenvalue. These curves cannot terminate since if they did, one 
could analytically continue from a region where the expression for / depends on a given 
dominant eigenvalue, to a region where it depends on a different eigenvalue. Owing to the 
Ug — i> 1/ug symmetry of the model, a given curve labelled by n approaches complex infinity 
in the same direction 6'„ as it approaches the origin. This yields (i) and (ii). From (H), it 
follows that 

< e„ < 7r/2 for < n < [s^ - 1/2] (10) 
which comprises [s^] values, and similarly, 

7r/2 < 0„ < TT for [s^ - 1/2] <n< [2s'^ - 1/2] (11) 

which again comprises [s^] values. Finally, if and only if s is half-integral, then the equation 
9n = 71 has a solution (for integral n), viz., n = (2s + l)(2s — l)/2. From the Ug u* 
symmetry of the phase diagram, the curve corresponding to this solution must lie on the 
negative real axis for all r, not just r — and r oo. This yields (iii)-(iv). □. 

We note that from (|^), the angular size of each region near the origin (or infinity) is 

A»=^ (12) 

Also, from the Ug —>■ u* symmetry of the phase diagram, it follows in particular, that for 
each curve starting out from the origin at 6'„, there is a complex conjugate curve at —On- As 
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is true of any theory at its lower critical dimensionality, the model is singular at K = oo, i.e, 
Us = 0. If and only if s is integral, the ra'th curve has another, m'th curve which is related 
to it hj 6m = On + (whence m = n + 2s^), so that as one travels through the origin on the 
n'th curve, one emerges on the other side on the m'th curve. 

As one can see from Fig. 1, as the boundary curves approach the unit circle \us\ = 1, some 
of them twist in s-dependent ways. Interestingly, several of the points where they cross the 
unit circle coincide with points which we inferred to be likely multiple (=intersection) points 
of the boundary curves of the complex-temperature phase diagrams for the corresponding 
spin s 2D Ising model. For example, for the ID s = 1 model, the n = and n = 1 curves cross 
the unit circle at Mi = z and Ui = e^*'^/^, respectively, and hence the complex-conjugate curves 
(n = 3 and n = 2) cross this circle at Ui = —i and Ui = e~^*'^/^. The points Ui = e^^*'^/^ are 
precisely the values which we inferred for the intersection points of boundary curves in the 
second ( "northwest" =NW) and third (SW) quadrants of the 2D s = 1 phase diagram from 
our calculation of partition function zeros For the ID s = 3/2 case, the n = 0, 1, 2 and 3 
curves cross the unit circle at the respective points M3/2 = e*'^/^, e^*'^/^, i, and e^*'^/^, and so 
forth for the complex conjugate curves. Among these, the points ±z and e^^*'^/^ are points 
inferred as likely intersection points of phase boundaries in the NW and SW quadrants of 
the M3/2 plane for the 2D s = 3/2 model 0]. Similar correspondences hold for s = 2. These 
are intriguing results. 

The spin 1/2 Ising model is equivalent to the two-state Potts model. While determining 
the CT phase diagram of the ID spin s Ising model, it is thus also of interest to compare 
this with that of the ID g-state Potts model, defined by Zp = J^a^ e"^''^^ with Hp = 
— JpJ2<nn'>^cT„a^, whcrc cr„ G {l,...,q}. We define Kp = [3Jp and up = e"^'' . The 
eigenvalues of the transfer matrix are A = u^^ — 1 (q' — 1 times) and A' = u'p^ + (g — 1). 
Setting Up = re*^, the equation for |A| = |A'| is 

g((g-2)r + 2cose) = (13) 

For q = 2, the solution is the imaginary axis, Up = ±ir, equivalent to the negative real 
axis in the ui/2 variable. For q ^ 2 (and g 7^ ), the solution is cos0 = — (g — 2)r/2 for 
< r < 2/(g-2), i.e., 

-1 + e^"^ 

up = (14) 
q-2 

for < < 27r, viz., a circle centered at up = — l/(g — 2) with radius l/(g — 2). Thus, for 
all q > 2, the complex-temperature phase diagram is qualitatively the same, consisting of a 



^Recall that the equivalence of Hp for g = 2 to (0) for s = 1/2 entails Kp = 2K. 
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PM phase containing the positive real u axis, and an O phase inside the circle This 
is quite different from the ID spin s Ising phase diagram, for which the number of phases 
does depend on s. This difference can be traced to the fact that the structure of the transfer 
matrix is qualitatively the same for different q in the Potts model (all diagonal entries equal 
to and all off-diagonal entries equal to 1), whereas it depends on s in the Ising model. 

In addition to the correspondences already above, we have found a further very interesting 
connection between the features of the ID spin s Ising model and certain properties which 
we have observed from our calculations of partition function zeros for the 2D (square-lattice) 
Ising model with s = 1, 3/2, and 2 One of the interesting features which we found in 
the 2D case was a certain number of finite arcs which protrude into the FM phase. (By the 
Us — > symmetry, there are also corresponding arcs protruding into the AFM phase). 
By combining our results with low-temperature series analyses in Refs. 0, 0], we observed 
that the divergences of the spontaneous magnetization M (which also imply divergences in 
the susceptibility |10|) occur at the endpoints of these arcs. The arc endpoints are thus of 



considerable interest for the study of CT singularities. 

For the cases s = 1, 3/2, and 2, the numbers of arcs protruding into and terminating in 
the first and second (NE and NW) quadrants of the (complex-temperature extension of the) 
FM phase are given in terms of the numbers of boundary curves of the corresponding ID 
spin s model as 

Ne,NE,2D = N,^NE,1D ' 1 = [s'] - 1 (15) 

and 

Ne,NW,2D = Nc,NW,lD = [s^] (16) 

SO that the total number of arcs protruding into the FM (equivalently, AFM) phase is 4[s^] —2. 
Our derivations in the present paper thus provide the explanation underlying our conjecture 
|Jl[| that the total number of arcs protruding into the FM (equivalently AFM) phase for the 
Ising model on the square lattice with arbitrary spin s is 

Ne,FM = Ne,AFM = max{0, 4[s^] - 2} (17) 

where we have incorporated also the exact result that Ng^pM = Ne^AFM = for s = 1/2. 

One can qualitatively describe how the complex-temperature phase diagram changes as 
one goes from d = 1 to d = 2 as follows. For d > 1, the model is analytic in the vicinity 
of Ms = (equivalently, it has a low-temperature series expansion with finite radius of 
convergence), so one knows that all of the curves which emanate from = must move 

"'One may also formally consider other values of q, as in the 2D Potts model. In particular, for q — Q, all 
eigenvalues are identically equal, and the phase diagram is trivial. 
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away from the origin. Evidently, the inner branches of the n = curve and its complex 
conjugate (n = — 1) curve move to the right and join at the point {us)c = e'^"^'^'^ 
which constitutes the critical point between the FM and PM phases; given the Us —* 1/us 
symmetry, this means that the outer branches no longer run separately to complex infinity 
but join each other at the critical point l/{us)c separating the PM and AFM phases of 
the 2D model. This leaves [s^] — 1 curves in the first (NE) quadrant, which form finite 
arcs; the inner and outer branches of these arcs protrude into and terminate in the FM and 
AFM phases respectively (these branches and phases are interchanged by the Ug — > 
mapping). The situation in the complex conjugate fourth (SE) quadrant is the same, by the 
Us — > u* symmetry of the phase diagram. In the second (NW) quadrant, the [s^] curves in 
ID correspond to the [s^] arcs in 2D, which again have inner and outer branches terminating 
in the FM and AFM phases. The symmetry under complex conjugation implies the same 
for the third (SW) quadrant. 

In summary, we have found an intriguing connection between the complex-temperature 
phase diagrams, as calculated exactly for ID and inferred from partition function zeros for 
2D, of the Ising model with higher spin. Our derivations here explain the basis for the 
conjecture that we made in Ref. [||] on the total number of arc endpoints in the FM phase 
and hence complex-temperature divergences in the magnetization for the square-lattice Ising 
model with arbitrary spin. 

This research was supported in part by the NSF grant PHY-93-09888. 
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Figure Captions 



Fig. 1. Complex-temperature phase diagrams in the Ug plane for the ID Ising model 
with (a) s = 1 (with u = Ui) and (b) s = 3/2 (with u = ^3/2)- The phase boundaries are 
shown as the dark curves (including the line along the negative Ug axis for s = 3/2). The 
unit circle is drawn in lightly for reference. 
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